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We investigate nonequilibrium entanglement generation in a chain of harmonic oscillators with 
time-dependent linear coupling. We use optimal control theory to determine the coupling modula- 
tion that leads to maximum logarithmic negativity for a pair of opposite oscillators and show that 
it corresponds to a synchronization of the eigenmodes of the chain. We further analytically relate 
the maximum attainable entanglement to the irreversible work done to produce it, thus bridging 
nonequilibrium entanglement production and nonequilibrium thermodynamics. 
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Entanglement is a puzzling feature of quantum the- 
ory. Entangled states are routinely created and manip- 
ulated in dozens of laboratories around the world: Con- 
trolled multiparticle entanglement has for example been 
achieved for up to eight qubits in ion traps [J 0]. In 
addition, entangled pairs of photons have been prepared 
over large distances of several kilometers in optical fibers 
[3, H[ as well as in free space d, @ . Due to their unique 
properties, entangled states have lead to many applica- 
tions in quantum communication and quantum informa- 
tion theory that are not possible at the classical level, 
such as quantum cryptography 0] and teleportation 0. 
And yet, the true nature of entanglement remains elusive, 
especially in many-body systems [9(. 

For quantum many-body systems at thermal equilib- 
rium, macroscopic entanglement has recently been re- 
lated to thermodynamic quantities like the magnetic sus- 
ceptibility Qi| , the internal energy [IJ, [l^] , and the heat 
capacity jl3l|. Some of these theoretical considerations 
have been successfully confirmed in experiments with 
quantum spin systems [H, 0, [H, . In the above ex- 
amples, equilibrium thermodynamic functions can thus 
be regarded as entanglement witnesses that can be used 
to detect and characterize entangled states. The situa- 
tion, however, is much different away from equilibrium 
and the question of how to relate entanglement produc- 
tion to irreversible thermodynamics is largely open. 

Our purpose in this paper is to determine the max- 
imum attainable entanglement generated in a nonequi- 
librium transformation and express the optimal value in 
terms of the irreversible work done during the process. 
The concept of irreversible work plays a fundamental role 
in the st udy of nonequilibrium thermodynamics of small 
systems [IjJ. To be specific, we shall consider a chain of 
quantum harmonic oscillators with time-dependent lin- 
ear coupling. This system is of direct relevance to the 
investigation of arrays of nanomechanical resonators [HI 
and ion crystals in multitrap arrangements [l9j . A no- 
table feature of the linear chain is that its dynamics is ex- 
actly solvable. Furthermore, since Gaussian states retain 
their Gaussian shape, we can rely on exact entanglement 
measures such as the logarithmic negativity [201 ] . 

The equilibrium entanglement properties of a closed 
harmonic chain with constant coupling have been studied 



in detail in Ref . (2l| . A remarkable result is that individ- 
ual oscillators are only entangled with their immediate 
neighbors. Contrary to classical correlations, entangle- 
ment hence does not occur over large distances within the 
chain. Clearly, this property also persists for adiabatic, 
that is, infinitely slow transformations of the Hamilto- 
nian, since the chain is then always in equilibrium. Re- 
cently, Eisert and coworkers have shown that a nonadi- 
abatic transformation - a sudden switch of the coupling 
strength - can significantly entangle even very distant 
oscillators, without having to control each of them indi- 
vidually [22| • They have considered the feasibility of such 
a scheme in the context of an array of double-clamped 
coupled nanomechanical oscillators and discussed the ex- 
perimental detection of the entanglement produced. 

In the following, we investigate nonequilibrium entan- 
glement generation in the harmonic chain from the point 
of view of irreversible thermodynamics. We use opti- 
mal control theory 23] to determine the modulation of 
the coupling coefficient that leads to maximum entangle- 
ment between distant pairs of oscillators. We find that 
optimal entanglement corresponds to a synchronization 
of the eigenfrequencies of the linear chain. Moreover, we 
show that the optimal value can greatly surpass that ob- 
tained from a sudden switch of the coupling strength. By 
further using recent findings on the nonequilibrium work 
distribution of a single oscillator [24|, we are able to re- 
late the optimal entanglement to the irreversible work 
performed during the nonadiabatic transformation. 

Harmonic chain. We consider a closed chain of N har- 
monic oscillators with identical frequency ujq and linear 
nearest-neighbour coupling with variable intensity c(i): 

1 N 

71=1 

To analyze the Hamiltonian (jTJ) , it is convenient to intro- 
duce the normal mode coordinates Q s and P s via, 
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The Hamiltonian of the chain can then be rewritten as 
a sum of independent oscillators with time-modulated 
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frequencies, w 2 (t) = w 2 + 4c(i) sin 2 (7rs/iV), 
i w 

where we have used the property, Q_ s = QJ, P_ s = PJ. 
The linear Heisenberg equations of motion for P s (t) and 
Q s {t) can readily be written down and solved [22 |. 

The Gaussian states of the linear chain are fully char- 
acterized by the symmetric 2N x 2N covariance matrix 
r, whose elements are defined by T qnqm = 2Ke(q n q m ), 
r ?„ Pm = 2Re(g n p m ) and T PnPm = 2Re(p n p m ). We 
have omitted the first moments, since they do not in- 
fluence the entanglement properties. Assuming for the 
time being that the chain is initially in the ground state, 
and that therefore T Qnqm (0) luq — T PnPm (0)/u>o — 6 n , m , 
T qnPm (0) = 0, the matrix elements can be written in 
terms of the second moments of the normal modes as, 

1 N 

(q n q m ) =JfYl e 2ms(n - m)/N (QsQl) , (4) 

s=l 

using (QrQl) = (QsQi)3r,s- Similar expressions are ob- 
tained for (q n p m ) and (p n p m )- The covariance matrix is 
hence entirely determined by the second moments of the 
normal coordinates. It is known that a time-dependent 
oscillator is squeezed when its frequency is changed nona- 
diabatically, while an adiabatic transformation leads to 
no squeezing (2f|. We can therefore express the quadra- 
tures of the normal oscillators in terms of the squeezing 
parameter r s and the rotation angle 6 S in the form, 

(QsQt) = 7T^(e" 2rs cos 2 9 S + e 2r ° sin 2 8 a ) , (5a) 
2lu s 

(P.pt) - ^(e^cos 2 ^ + e ^sin 2 s ) , (5b) 

(QsPl) = sinh(2r s )sin0 s cos0 s . (5c) 

The time dependence of u) s , r s and 9 S is controlled by 
the linear coupling coefficient c(t). Standard thermo- 
dynamics considers infinitesimal changes that take place 
infinitely slowly. In the spirit of finite-time thermody- 
namics [26j. we will in the sequel consider infinitesimal 
changes that can occur at arbitrary rate. We will thus as- 
sume that the coupling intensity is small, c(t) <C Wq, but 
otherwise arbitrary. Moreover, without loss of generality, 
we will set the frequency wq = 1. 

Maximum entanglement. The entanglement between 
any two oscillators n and m in the chain can be charac- 
terized with the help of the logarithmic negativity [2(| , 

i 4 

£W = --5>g 2 [Min(l,|^|)] . (6) 

»=i 

Here, Zj are the symplectic eigenvalues of the covariance 
matrix of the partially transposed reduced density oper- 
ator of the two oscillators. In the present investigation, 



we choose two opposite oscillators, m = N/2 + n. The 
latter configuration corresponds to the largest possible 
distance in the chain. It can also be shown to corre- 
spond to the largest value of the logarithmic negativ- 
ity due to symmetry reasons. It is worth noting that 
in this case the exponential factor in Eq. (QJ simplifies 
to e 2ms{n-m)/N _ (— l) s 5 thus distinguishing even and 
odd normal modes. Furthermore, in view of the transla- 
tional invariance along the chain, the logarithmic nega- 
tivity does not depend on the particular value of n. 

The symplectic eigenvalues of the 4x4 reduced co- 
variance matrix T m ^ n can be computed explicitly using 
the parametrization ([5]). A careful examination of their 
structure reveals that for moderate squeezing the loga- 
rithmic negativity ^ is maximum when the angles 6 S 
take special values. These optimal angles are given by 
9 S = 7r/4(mod7r) for s even and 9 S = 37r/4(mod7r) for 
s odd. Hence, for fixed squeezing, maximum entangle- 
ment is achieved when both even and odd oscillators 
are synchronized. From Eqs. (JSJ), it follows that, at 
the optimal angles, (Q s Qt) = cosh(2r s )/2w s , (P s Ps) — 
cosh(2r s )w s /2 and (Q s Pj) = (~l) s sinh(2r s )/2. Com- 
bining Eqs. ^ and ([5]), with lu s ~ uo, we then find that 

(in) = iPn) = a and (QnQm) = (PnPm) = b. As a result, 

the symplectic eigenvalues in Eq. ^ are simply given by, 

\h\ = y/[aTb-(cTd)][a±b + c±d\ , (7) 

with c = (q n Pn) and d — (q n Pm)- We therefore obtain 
that the maximum attainable logarithmic negativity is: 




Equation ([8]) shows that entanglement is zero for van- 
ishing squeezing and increases with increasing squeez- 
ing. By comparing the simplified symplectic eigenval- 
ues ([7]) with their general expression, we further find 
that the optimal angles are valid provided the condition 
cJ2 s exp(2r s )/4 <C N is satisfied. As we shall discuss 
below, this condition is not very restrictive. 

Optimal control theory. Optimal control theory (OCT) 
is a powerful method that allows to determine the func- 
tion that optimizes a given cost functional [23|. It 
can be regarded as a generalization of the calculus 
of variations and has a long tradition in finite-time 
thermodynamics (26j. We apply Pontryagin's princi- 
ple to find the function c(t) that optimizes the entan- 
glement between opposite oscillators. To this end, we 
take as the cost functional to minimize, the argument 
of the maximum attainable logarithmic negativity |8|). 
Es, ? dd ex P(- 2r s)E s , e j, e „exp(-2r s ). It is clear that the 
minimum of this quantity will lead to maximum entangle- 
ment. The constraints given by the equations of motion 
of the normal modes ([2|) enter the optimization prob- 
lem in the form of Lagrange multipliers. The optimal 
function c(t) is eventually obtained by steepest descent. 
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FIG. 1: Logarithmic negativity, Eq. (B). for opposite oscil- 
lators as a function of time for different modulations of the 
coupling strength c(t) between c = and c = 0.05. The lower 
shaded curve corresponds to a sudden switch at t = 0, while 
the curve in the center is generated by the optimal modula- 
tion shown in Fig. [3] The two upper curves are the maximum 
attainable entanglement for T = and T = 0.5. 



The results of the numerical implementation of optimal 
control theory for a chain of N = 8 oscillators are sum- 
marized in Figs. [U [2] and [3] 

Figure [I] shows the logarithmic negativity Em of oscil- 
lators 1 and 5 as a function of time for various coupling 
intensities c(t). The shaded curve at the bottom corre- 
sponds to a sudden switch from c = to c = 0.05. The 
oscillatory curve in the center depicts the entanglement 
generated using OCT for the same parameters. The up- 
per curve, on the other hand, is the maximum attainable 
entanglement at zero temperature, Eq. (J8J), when the an- 
gles 8 S of the normal oscillators are perfectly synchro- 
nized at all times. We observe that the entanglement 
produced by the sudden switch is one order of magni- 
tude smaller than the maximum attainable value in the 
present situation. Moreover, we note that the optimal 
logarithmic negativity (almost) reaches the maximum at- 
tainable value E' N iax . Figure [5] shows the angles of the 
normal oscillators, corresponding to the optimal modula- 
tion, in two short time intervals: the first around the peak 
of highest entanglement (denoted by B in Fig. [I]) and the 
second around a minimum of the logarithmic negativity 
(interval A). The (nearly perfect) synchronization of the 
even and odd modes at the center of the interval is clearly 
visible in case B, while it is absent in case A. The optimal 
modulation c{t) is plotted in Fig. [3] and consists of three 
different parts. The first part is a sequence of square 
pulses between c — and c = 0.05. The latter is the re- 
sult of the optimization procedure that tries to maximally 
squeeze each normal oscillator in order to minimize the 
argument of Eq. (JSJ) ; the maximal entanglement E^ ax is 
determined by the duration of this sequence. In the sec- 
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FIG. 2: Angles 9 3 of the normal modes as a function of time 
for the two intervals A and B shown in Fig. [T] (due to the 
symmetry of the chain, angles are pairwise equal). The syn- 
chronization of even and odd modes is clearly visible at the 
center of interval B and completely absent in interval A. 

ond part, the coupling c(t) is kept constant allowing the 
angles to synchronize [27| ; the synchronization time is di- 
rectly controlled by the value of c: a larger coupling will 
lead to a faster synchronization. Since each normal mode 
evolves with a different frequency, the synchronization of 
the angles cannot persist in time. However, by stopping 
the modulation of c(t) (dashed line in Fig. [3]), once the 
biggest entanglement has been attained, the peak value 
can be maintained (see dashed line in Fig. [IJ . 

Dissipative work. We next relate the nonequilibrium 
entanglement produced by a nonadiabatic switching of 
the coupling to the average dissipated work. The dissi- 
pated work is a central concept in thermodynamics and 
is intim ately connected to the irreversible entropy pro- 
duction [13|; it is defined as the difference between the 
total work and the reversible work (or free energy differ- 
ence), Wdis = W — AF. For an adiabatic transformation, 
Wdis = and the total work is given by AF. The dissi- 
pated work for a frequency-modulated normal oscillator 
of the chain can be simply expressed in terms of the fre- 
quency and the squeezing parameter at time t [28(: 

W d is,s = w s sinh 2 r s . (9) 

The above equation indicates that there is a one-to-one 
correspondence between the dissipated work Wdis.s and 
the squeezing parameter r s : small squeezing yields small 
dissipated work, Wdis.s — u s r1, whereas large squeez- 
ing leads to an exponential increase of the dissipated 
work, Wdis,s — &s exp(2r s )/4. The total work dissipated 
in the chain is given by the sum over all eigenmodes, 
Wdis = J2s W dis, s - By combining Eqs. © and ©, we 
can now directly relate the maximum attainable loga- 
rithmic negativity to the irreversible work required to 
generate it. This relationship is of importance for sev- 
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FIG. 3: Optimal coupling modulation c(t) that generates the 
logarithmic negativity shown in the center of Fig. [1] 

eral reasons: First, from a fundamental point of view, 
it bridges nonequilibrium entanglement production with 
nonequilibrium thermodynamics; second, it gives us the 
maximum amount of entanglement that can be produced 
by dissipating a given amount of work, in the ideal situ- 
ation where all the phases are synchronized. One could 
for instance dissipate considerably more work by highly 
squeezing few oscillators in the chain, irrespective of their 
angles, and produce entanglement far below the maxi- 
mum possible value. In that respect, this relationship is 
similar in spirit to the Carnot efficiency of a heat engine 
that gives the maximum efficiency that can be achieved 
for given temperatures, in an ideal cycle. 

Effect of temperature. Our discussion so far was limited 
to a chain at zero temperature. Finite temperatures can 



easily be taken into account by multiplying the initial co- 
variance matrix r(0) by the factor f(T) = coth(w /2T). 
The new symplectic eigenvalues are then merely multi- 
plied by the factor f(T). As a consequence, the max- 
imum attainable logarithmic negativity is still given by 
its expression (|8[) . with the argument now multiplied by 
f{T). Since f(T) > 1, temperature reduces the maximal 
value of entanglement, as expected. Figure [1] shows the 
maximal logarithmic negativity for T = 0.5. A good es- 
timate of the highest temperature T m at which entangle- 
ment vanishes can be obtained by replacing the squeezing 
parameters of the different eigenmodes by the average 
squeezing, R = J2 S r s/N. We then find that entangle- 
ment survives up to T m = Wo/ln(coth(2iJ)). 

Discussion. One of the main limitations that have to 
be considered in real systems is the validity of the har- 
monic approximation. As the oscillators are squeezed, 
their energy increases and anharmonicities may appear. 
For the optimal modulation of Fig. [3] with a squeezing 
sequence up to t = 20, R = 0.4 and Wdis = 1-2 (in units 
of h) . This leads to an increase of 30% of the total energy 
of the chain from E = 4 to E + Wdis = 5.2. A tenfold 
energy increase is achieved with a squeezing sequence up 
to t = 85, yielding R = 1 and E^ ax = 1.6. These param- 
eters are well within the bounds of the criterion given 
below Eq. ([5]). A second limitation is environment in- 
duced decoherence. Fortunately, as shown in Ref. [2^ |. 
entanglement in the chain is particularly robust against 
external noise sources at low temperatures. 
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